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Abstract 

The Dalitz plot distribution of B° —> D°K + ir~ decays is studied using a data sample 
corresponding to 3.0 fb -1 of pp collision data recorded by the LHCb experiment 
during 2011 and 2012. The data are described by an amplitude model that contains 
contributions from intermediate A'*(892)°, A'*(1410)°, A'2(1430)° and -D?;(2460)~ 
resonances. The model also contains components to describe broad structures, 
including the A"q( 1430)° and Aq( 2400)~ resonances, in the Kir S-wave and the Dir S- 
and P-waves. The masses and widths of the A>q( 2400)“ and A>2(2460) _ resonances 
are measured, as are the complex amplitudes and fit fractions for all components 
included in the amplitude model. The model obtained will be an integral part 
of a future determination of the angle 7 of the CKM quark mixing matrix using 
B° —» DK + ir~ decays. 
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1 Introduction 


Dalitz plot (DP) analysis of B° —y DK + n~ decays has been proposed as a way to measure 
the unitarity triangle angle 7 |l,2|. The sensitivity to 7 = arg [— V ud V* b /(V cd V* b )\, where 
V xy are elements of the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix |3,[Ij, 
originates from the interference of b — y cus and b —y ucs amplitudes. Such interference 
occurs when the neutral D meson is reconstructed in a final state that is accessible to both 
D° and D° decays and therefore corresponds to an admixture of the two states [ 5 ,61. One 
of the largest components of the B° —* DK + tt~ final state is B° —> DK* (892)°, for which 
both b —* cus and b —* ucs amplitudes are colour suppressed, making them comparable 
in magnitude and potentially enhancing CP violation effects |7|. Decay diagrams for the 
quasi-two-body contributions from B° —>■ DK*{ 892)° and B° —> K + decays 

are shown in Fig. [l] Observables sensitive to 7 have been measured by LHCb with a 
quasi-two-body approach |&|, but a DP analysis is expected to be more sensitive because 
interference between resonances provides the possibility to resolve ambiguities in the 
determination of 7. 

In order to determine 7 with this method, it is necessary to have an amplitude model 
of the B° —y D°K + ti~ decay that proceeds through the favoured b —» cus transition. This 
can be achieved by reconstructing the D° meson through the K + n~ decay; in this way 
the contribution from the b —» cus amplitude is dominant and effects due to b —> ucs 
amplitudes can be neglected. While the analysis of this decay chain is not itself sensitive to 
7, its outcome will be an integral part of a future analysis using, for example, D —s- K + K~ 
decays where CP violation effects are expected as the final state is common to both D° 
and D° decays. 

Dalitz plot analyses of B meson decays to final states containing a charmed meson and 
two charged particles (either pions or kaons) also provide opportunities for studies of the 
spectroscopy of charmed mesons. Results in this area have recently been obtained from DP 
analyses of B° s —>• D°K~ 7r + [9,10 , B + —y D~K + n + 111 and B° —> 12 decays, 

all from LHCb. As the branching fraction for B° —> D°K + tt~ decays is smaller than that 
for B° —> D 0 7i + ir~ decays, the analysis presented in this paper is not as sensitive to the 
parameters of charm resonances as that of Ref. jl2|. However, the much larger sample of B 



Figure 1: Decay diagrams for the quasi-two-body contributions to B° —> DK + ir from (a) 
B° -> DK* (892)° and (b) B° -> Z^(2460)-A+ decays. 
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mesons available at LHCb compared to that used in the only other published DP analysis 
of B° —y D°7 i + t:~ decays from the Belle collaboration 13 allows useful results on excited 


charm mesons to be obtained. Moreover, results on charm meson spectroscopy obtained 
from DP analysis of B° —>• D°K + n~ decays provide important independent cross-checks 
of results from studies of the B° — y D 0 7i + n~ DP, as possible biases due to other structures 
in the Dalitz plots are different between the two modes. 

In this paper an amplitude analysis of the B° —> D°K + tt~ decay is reported. The 
inclusion of charge-conjugate processes is implied throughout the paper. The analysis 
is based on a data sample corresponding to an integrated luminosity of 3.0 ftV 1 of pp 
collision data collected with the LHCb detector in 2011 and 2012, when the collision 
centre-of-mass energy was ^fs = 7TeV (l.Ofb -1 ) and yfs = 8TeV (2.0fb -1 ), respectively. 


Previously, the branching fraction for the three-body decay has been measured 14,15 
and the B° —> D°K*(892)° 16,17 and B° —>■ D\ (2460) ~K + 14 contributions have been 


measured using quasi-two-body approaches: however this is the first DP analysis of the 
B° D°K+7 T" decay. 

The paper is organised as follows. A description of the LHCb detector, reconstruction 
and simulation software is given in Sec. [2} The selection of signal candidates is described 
in Sec. [3], and the determination of signal and background yields is presented in Sec. |4j 
An overview of the Dalitz plot analysis formalism is given in Sec. [5j and details of the 
implementation of the amplitude analysis are presented in Sec. [6| The evaluation of 
systematic uncertainties is described in Sec. [7], with results and a brief summary given in 
Sec. [U 


2 LHCb detector and software 


The LHCb detector |18jjl9| is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < rj < 5, designed for the study of particles containing b or c quarks. The detector 
includes a high-precision tracking system consisting of a silicon-strip vertex detector |20 
surrounding the pp interaction region, a large-area silicon-strip detector located upstream 
of a dipole magnet with a bending power of about 4 Tin, and three stations of silicon-strip 
detectors and straw drift tubes [21 placed downstream of the magnet. The polarity of 


the dipole magnet is reversed periodically throughout data-taking. The tracking system 
provides a measurement of momentum, p, of charged particles with a relative uncertainty 
that varies from 0.5% at low momentum to 1.0% at 200 GeV (natural units with c — h—1 
are used). The minimum distance of a track to a primary vertex, the impact parameter 
(IP), is measured with a resolution of (15 + 29/px) pm, where pr is the component of 
the momentum transverse to the beam, in GeV. Different types of charged hadrons 
are distinguished using information from two ring-imaging Cherenkov detectors [221. 
Photon, electron and hadron candidates are identified by a calorimeter system consisting of 
scintillating-pad and preshower detectors, an electromagnetic calorimeter and a hadronic 
calorimeter. Muons are identified by a system composed of alternating layers of iron and 
multiwire proportional chambers 


23 
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The trigger 124 consists of a hardware stage, based on information from the calorimeter 
and muon systems, followed by a software stage, in which all tracks with p t > 500 (300) MeV 
are reconstructed for data collected in 2011 (2012). The software trigger used in this 
analysis requires a two-, three- or four-track secondary vertex with significant displacement 
from any primary pp interaction vertex (PV). At least one charged particle must have 
Pt > 1.7 GeV and be inconsistent with originating from a PV. A multivariate algorithm 25 


is used for the identification of secondary vertices consistent with the decay of a b hadron. 

In the offline selection, the objects that prompted a positive trigger decision are 
associated with reconstructed particles. Selection requirements can therefore be made not 
only on whether the hardware trigger decision was due to a signature in the calorimeters 
or in the muon system, but on whether the decision was due to the signal candidate, other 
particles produced in the pp collision, or a combination of both. Signal candidates are 
accepted offline if at least one of the final-state particles created a cluster in the hadronic 
calorimeter with sufficient transverse energy to fire the hardware trigger. Events that are 
triggered at the hardware level by another particle in the event are also retained. After 
all selection requirements are imposed, 67 % of events in the sample were triggered by 
the decay products of the signal candidate, while the remainder were triggered only by 
another particle in the event. 

Simulated events are used to characterise the detector response to signal and certain 


types of background events. In the simulation, pp collisions are generated using Pythia 26 
with a specific LHCb configuration 


27 


by EvtGen 28 , in which final-state radiation is generated using PHOTOS 29 


Decays of hadronic particles are described 

The 


interaction of the generated particles with the detector, and its response, are implemented 
using the Geant 4 toolkit 


30 as described in Ref. 31 


3 Selection requirements 


The selection requirements follow closely those used in Ref. 
B° -A 5° 7T 
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7T 


The more copious 

decay is topologically and kinematically similar to the B° —>■ D°K + ti~ 
channel, allowing it to be used as a control mode to optimise the selection requirements. 
Loose initial requirements are used to obtain a visible signal peak of D q ti + ti^ candidates. 
The tracks are required to be of good quality and must be above thresholds in p, p T and 
Xj P , where Xip is defined as the difference in y 2 of a given PV reconstructed with and 
without the considered particle. The D° —> K + n~ candidate must satisfy criteria on its 
vertex quality (y 2 tx ) an d flight distance from any PV and from the B candidate vertex, 
and must have invariant mass m(K + 7T~ ) in the range 1814-1914 MeV. A requirement- 
on the output of a boosted decision tree that identifies D° —> K + n~ decays originating 
from b hadron decays (D° BDT) 


32 


33 


is also applied. Candidate B mesons are selected 
with requirements on invariant mass, Xip and on the cosine of the angle between the B 
momentum vector and the line from the PV under consideration to the B vertex (cos#dir)- 
A requirement is placed on the y 2 of a kinematic fit 34 to the candidate’s decay chain in 
which the D° mass is constrained to its nominal value. The four final-state tracks must 
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satisfy pion and kaon particle identification (PID) criteria. 

A neural network 35 is used to discriminate between signal decays and combinatorial 
background. The sPlot technique [36 J, with the B candidate mass as the discriminating 
variable, is used to separate statistically B° —» D 0 tt + tt~ decays from background. The 
signal and background weights returned by this method are applied to the candidates, 
which are then used to train the network. The network is trained using 16 variables. They 
include the \ip of the four final-state tracks and the following variables associated to 
the D° candidate: yjp, y^ tx , the square of the flight distance from the PV divided by 
its uncertainty squared (x| ight ), cos^dir, and the output of the D° BDT. In addition, the 
following variables associated to the B candidate are included: pt, Xjp, Xvtx; Xm g ht> an d 
cos^dir- The px asymmetry 10 and track multiplicity in a cone with half-angle of 1.5 


units of the plane of pseudorapidity and azimuthal angle (measured in radians) around 
the B candidate flight direction, which contain information about the isolation of the B 
candidate from the rest of the event, are also used. The input quantities to the neural 
network depend only weakly on the position of the candidate in the B decay Dalitz plot 
and therefore any requirement on the network output cannot appreciably bias the DP 
distribution. A requirement imposed on the network output reduces the combinatorial 
background remaining after the initial selection by a factor of five while retaining more 
than 90 % of the signal. 

The B° —> D 0 K + n~ candidates must satisfy the same selection as the B° —> D 0 7i + 7r~ 
decays, except for the PID requirement on the positively charged track from the B decay 
vertex, which is imposed to preferentially select kaons rather than pions. The PID efficiency 
is evaluated from calibration samples of D° —> K~n + decays from the D* + —> D° ti + decay 
chain. The kinematics of this decay chain can be exploited to obtain clean samples without 

The PID efficiency of the requirements on the four tracks 
3 and varies depending on the kinematics of the tracks, as 


using the PID information 22 


in the final state is around 50 / 
described in detail in Sec. 16.11 

Candidates are vetoed when the difference between m(K + n~n~) and m{K + n~) lies 
within ±2.5 MeV of the known D*(2010) _ -D° mass difference 37 to remove background 
containing D*~ — > D°tt~ decays. Candidates are also rejected if a similar mass difference 
calculated with the pion mass hypothesis applied to the bachelor kaon satisfies the same 
criterion. To reject backgrounds from B° — > D~K + , D~ —> K + 'k~ti~ decays, it is required 
that the combination of the pion from the D° candidate together with the two bachelor 
particles does not have an invariant mass in the range 1850-1890 MeV. Additionally, 
candidates are removed if the pion and kaon originating directly from the B° decay 
combine to give an invariant mass consistent with that of the D° meson (1835-1880 MeV). 
This removes candidates with the D° wrongly reconstructed as well as potential background 
from B° —>- D°D° decays. Other incorrectly reconstructed candidates are removed by 
vetoing candidates where the pion from the D° decay and the kaon originating directly 
from the B° decay give an invariant mass in the range 1850-1885 MeV. At least one of 
the pion candidates is required to have no associated hits in the muon system to remove 
potential background from B° — > J/^K* 0 decays. Charmless decays of b hadrons are 
suppressed by the use of the D° BDT and further reduced to a negligible level by requiring 
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that the D° candidate vertex is separated from the B° decay vertex by at least 1 mm. 

Signal candidates are retained for further analysis if they have invariant mass in the 
range 5100-5900 MeV. After all selection requirements are applied, fewer than 1 % of 
selected events also contain a second candidate. Such multiple candidates are retained and 
treated in the same manner as other candidates; the associated systematic uncertainty is 
negligible. 


4 Determination of signal and background yields 


The signal and background yields are determined from an extended maximum likelihood 
fit to the B candidate invariant mass distribution. The fit allows for signal decays, 
combinatorial background and contributions from other b hadron decays. The decay chain 
B° —> D*°K + n~ with D*° —> D °7 or D 0 7i° forms a partially reconstructed background 
that peaks at low B candidate mass as the neutral particle is not included in the candidate. 
Misidentified B° —> D^°n + TT~, A® —>• D^°K + p and Djh —)■ D^°K + K~ decays are found 
to contribute to the background. A possible contribution from the highly suppressed 
B° s —>• D°K + ir~ decay mode is also included in the fit. 

The signal peak is modelled with the sum of two Crystal Ball 


38 functions, which 


have tails on opposite sides and which share a common mean. The tail parameters are 
fixed to the values found in fits to simulated signal decays. The relative normalisation 
and the ratio of widths of the two functions are constrained, within uncertainties, to 
the values determined in a fit to B° —s- D° 7 r + 7 r - data. The B° s — D°K + tt~ shape is 
modelled identically to the signal peak, with the difference between the B° and B° s mass 

An exponential shape is used to model the 


37 


peak positions fixed to its known value 
combinatorial background. 

Misidentified and partially reconstructed backgrounds are modelled with smoothed 
non-parametric functions. Simulated samples are used to obtain the shape for B° —> 
D* 0 K + n~ decays, with the D*° —> D °7 and D 0/ n 0 contributions generated in the correct 
proportions |37|. To account for differences between simulation and data, for example 
in the polarisation of the D*° meson, the B candidate invariant mass distribution for 
partially reconstructed B° —» D*°K + n~ decays is allowed to be shifted by an offset that 
is a free parameter of the fit. The misidentified background shapes are also obtained 

decays the D° and D*° contributions are 
The backgrounds from 


7r 1 7r 


37 


from simulated samples. For B° —>■ ZA*)° 
combined according to their relative branching fractions 

A® —> D^°K + p and —> D^°K + K~ decays are assumed to have equal branching 

fractions for the decays with D° and D*° mesons in the final state, as the decays involving 
D*° mesons have not yet been measured. The simulated samples are reweighted according 
to the relevant particle identification and misidentification probabilities and to match known 
DP distributions 12,39-41 . The (mis)identification probabilities take account of track 
kinematics and are calculated using calibration samples of D* + —> D°ti + , D° — > K~7i + 
and A —>■ pir~ decays 19,22 . The yields of the misidentified backgrounds relative to the 


signal are constrained, within uncertainties, to their expected values based on their known 
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Table 1: Yields from the fit to the D°K + ir data sample. The full mass range is 5100-5900 MeV 
and the signal region is 5248.55-5309.05 MeV. 


Component 

Full mass range 

Signal region 

B° ->■ D°K + 7t- 

2576 ± 72 

2344 ± 66 

B° s D°K + 7T" 

55 ±27 

1 ± 1 

comb. bkgd. 

5540 ± 187 

684 ± 23 

B° D*°K + tt- 

1750 ± 99 

6 ± 1 

B° _+ B (*)o n + 7T - 

485 ± 47 

51 ±5 

jo £(*)o K+p 

95 ±26 

18 ±5 

B o d(*)or+K~ 

127 ±27 

10 ±2 

B o d(*)°k+K~ 

54 ± 18 

14 ±5 




Figure 2: Results of the fit to the B candidate invariant mass distribution with (a) linear and 
(b) logarithmic y-axis scales. The components are as described in the legend. 


branching fractions [12,39,401 and m ^ identification probabilities. 

There are 14 free parameters in the fit model: the mean and width of the signal shape, 
the relative normalisation and relative width of the two Crystal Ball functions, the slope 
of the exponential function, the offset of the B° —y D*°K + tt~ shape and the yields of the 
eight contributions. The result of the fit is shown in Fig. [2] and the yields are summarised 
in Table [lj The yields are also reported within the signal region used in the Dalitz plot- 
fit, corresponding to ±2.5 widths of the B° signal shape (5248.55-5309.05MeV). The 
distribution of candidates in the signal region over the Dalitz plot is shown in Fig. !»)• 
Note that a B° mass constraint is applied to calculate the variables that are used to 
describe the Dalitz plot 42 , improving the resolution of those variables and giving a 


unique kinematic boundary. 


6 



















12 

> 

0) 

o 10 


+ 


1 1 ' ' 1 ' 1 1 1 ' 1 1 ' ' 1 1 1 1 1 ' 

CL 1 

: 1 — 1 — 1 —t !— 1 —i— 1 —'— 1 —i— 1 —'— 1 —i— 1 —'— 1 : 

• LHCb 

0.9 

i- . .• _ ._____ -= 

4 . (a) 

0.8 

“ — 


0.7 

— • — 

■ 

0.6 



0.5 

h T&T -j 

: T.#V . ; 

0.4 

!- -i 


0.3 

I- 4 

- : ^ 

0.2 

L LHCb J 

- t;.. . _ . 

0.1 

7 (b) 7 

_1_1_._1_1_._1_^_,_._I_1_._._1_1_1_1_1_._ 

0 

_._,_l_^^_i_^._i_■_._,_i_._._^ 


10 15 20 

m 2 (D\-) [GeV 2 ] 


0.2 


0.4 


0.6 


0.8 


Figure 3: Distribution of B° —> D°K + ir candidates in the signal region over (a) the Dalitz plot 


and (b) the square Dalitz plot. The definition of the square Dalitz plot is given in Sec. 6.1 


5 Dalitz plot analysis formalism 


In B° — D°K + decays, resonances are expected in the m 2 (D°7i ~) and m 2 {K + Ti ~) 
combinations, so the Dalitz plot shown in Fig. [3]^a) is defined in terms of these two 
invariant mass squared terms. For a fixed B° mass, these two invariant mass squared 
combinations can be used to calculate all other relevant kinematic quantities. 

The isobar model 43-45 is used to describe the complex decay amplitude. The total 


amplitude is given by the coherent sum of amplitudes from resonant and nonresonant 
intermediate contributions and is given by 


N 

A (m 2 (-D°7r~), m 2 (K + n~)) — CjFj (m 2 (D°7r~), m 2 (K + n~)) , (1) 

j = 1 

where c 3 are complex coefficients describing the relative contribution for each intermediate 
process. The resonant dynamics are contained in the Fj (m 2 (D 0 TT~),m 2 (K + 7T~)^ terms 
that are normalised such that the integral of the squared magnitude over the DP is unity 
for each term. For a D°ti~ resonance Fj (^m 2 , m 2 ( y K + ir~)^ is given by 

F (m 2 (D°7 r“),m 2 (/t + 7 r“)) = R (m(D°n~)) x X(\p\ r BW ) xX(|g|r BW ) xT(p,q), (2) 


where p is the bachelor particle momentum and q is the momentum of one of the resonance 
daughters, both evaluated in the D 0 ti~ rest frame. The functions R, X and T are described 
below. 

The X(z) terms are Blatt-Weisskopf barrier factors [46], where z = \q\ r B w or \p\ r BW 
and r B -w is the barrier radius which is set to 4.0 GeV -1 ~ 0.8 fm 10 for all resonances. 
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The barrier factors are angular momentum dependent and are given by 


L = 0 : X(z) 
L = 1 : X(z) 
L = 2 : X(z) 
L = 3 : X(z) 


l+z 2 


Zq-\-3zq-\-9 

z 4 +3z 2 +9 


Zg + 62 q +45zq +225 
z 6 +6z 4 +45z 2 +225 


( 3 ) 


where Zq is the value of z at the pole mass of the resonance and L is the orbital angular 
momentum between the resonance and the bachelor particle. Since the parent and daughter 
particles all have zero spin, L is also the spin of the resonance. 


The T(p, q) terms describe the angular distributions in the Zemach tensor formalism 47 


48 and are given by 


L = 0 
L = 1 
L = 2 
L = 3 


T(p,q) = 1, 

T(p,q) = -2 p-q, 

T(p,q) = § [3(p- q) 2 - (\p\\q\) 2 } , 

T(p,q) = ~ ff [5(p- q) 3 - 3{p ■ q)(\p\\q\) 2 ] 


(4) 


These expressions are proportional to the Legendre polynomials, Pl(x), where x is the 
cosine of the angle between p and q (referred to as the hclicity angle). 

The R (m(-D°7r - )) functions are the mass lineshapes. Resonant contributions are 
typically described by the relativistic Breit-Wigner (RBW) function 


R(m ) = 


(mg — m 2 ) — % moT(m) 
where the mass-dependent decay width is 


(5) 


T(m) = T 0 M- 


(:) 


2L+1 


(—)x 2 (qr BW ), 

\ m / 


( 6 ) 


where qo is the value of q — \q \ when the invariant mass is equal to the pole mass of the 
resonance, mo. 

The large phase space available in B decays allows for the presence of nonresonant 
amplitudes (i.e. contributions that do not proceed via a known resonance) that vary across 
the Dalitz plot. An exponential form factor (EFF) has been found to describe nonresonant 


contributions well in several DP analyses of B decays 49 


R{m) = exp [—am 2 ] 


(7) 


where a is a shape parameter that must be determined from the data and m is a two-body 
invariant mass (m(D 0 7i~) in this example). 

The RBW function is a good model for narrow resonances that are well separated from 
any other resonant or nonresonant contribution of the same spin. This approach is known 
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to break down in the Ktc S-wave because the Jig(1430) resonance interferes strongly with 
a slowly varying nonresonant term, as described in Ref. j50 . The LASS lineshape 51 has 
been developed to combine these two contributions, 


R{m) = 


m 


q cot 5b — iq 


exp [2 i5 


B 


™ 0 r 0 ^ 


{ml -in?) -im 0 r 0 ±^ 


, f 1 1 

where coto# =-1— rq, 

aq 2 


( 8 ) 

(9) 


and where mo and To are the pole mass and width of the /Lq( 1430) state, and a and r are 
shape parameters. 

The Dn S-wave nonresonant contribution can be described by the “dabba” line- 
shape 52 , defined as 

R{m ) = 

where 


B'(m 2 )(m 2 — sa)p 


1 — f3(m 2 — mhj — iB'(m 2 )(m 2 — sa)p ’ 


( 10 ) 


B'{m 2 ) = b exp [— oi{m 2 — m^ in )] . (11) 

Here m m ; n is the invariant mass at threshold, = m 2 D — 0.5 m 2 is the Adler zero, p is 
a phase-space factor and b , a and f3 are parameters with values fixed to 24.49 GeV -2 , 


0.1 GeV 2 and 0.1 GeV 2 , respectively, according to Ref. 52 


T -2 


Ignoring reconstruction and selection effects, the DP probability density function would 


be 


■Pphys (m 2 (D°7T ), m 2 (A' + 7r )) = 


|^4(m 2 (D°7r ),m 2 (iP + 7r )) 


( 12 ) 


//dp l^l 2 dm 2 (D°7 r-)dm 2 (A'+7r- 

where the dependence of A on the DP position has been suppressed in the denominator for 
brevity. The primary results of most Dalitz plot analyses are the complex coefficients, as 
defined in Eq. ([!]). These, however, depend on the choice of phase convention, amplitude 
formalism and normalisation used in each analysis. The convention-independent quantities 
of fit fractions and interference fit fractions provide a way to reliably compare results 
between different analyses. Fit fractions are defined as the integral of the amplitude for 
each intermediate component squared, divided by that of the coherent matrix element 
squared for all intermediate contributions, 

12 


AA,= 


ffj 


DP 


\ C jFj 


m 2 (D°n ),m 2 (K + n ))|" dm 2 (D°n )dm 2 (K + n ) 


//dp 1*4 dm 2 (D°n~) dm 2 (K + n~) 


(13) 


The fit fractions need not sum to unity due to possible net constructive or destructive 
interference, described by interference fit fractions defined (for i < j only) by 


FF ■ — 

11 u 


ff DP 2 7Ze \ciCjFiF*] dm 2 (D°n )dm 2 (K + n 
// dp |*4| 2 dm 2 (D°TT~) dm 2 (K + TT~) 


( 14 ) 


where the dependence of F} ; and A on the DP position has been omitted. 
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Figure 4: Efficiency variation as a function of SDP position for candidates triggered by (a) signal 
decay products and (b) by the rest of the event. The vertical white stripe is due to the D* veto 
and the curved white band is due to the D° veto. 


6 Dalitz plot fit 

6.1 Signal efficiency 

The variation of the efficiency across the phase space is studied in terms of the square 
Dalitz plot (SDP). The SDP is defined by variables m! and 9' that range between 0 and 1 
and are given for the D°K + ti~ case by 

1 ( m(D°n~) — m"d> n _ \ l __ 

m! = — arccos [ 2--1 ] and 9' = —9(D°tt ), (15) 

n \ m max_ _ m mm _ I 7T 


where = mgo — mx+ and = m^ 0 + m^- are the kinematic boundaries of 

m(D°7i~) allowed in the B° —y D°K + n~ decay and 9(D°n~) is the helicity angle of the 
D 0 7i~ system (the angle between the K + and the D° meson in the D 0 7i~ rest frame). The 
distribution of selected events across the SDP is shown in Fig. |3](b). 

Efficiency variation across the SDP is caused by the detector acceptance and by 
trigger, selection and PID requirements. The efficiency is evaluated with simulated samples 
generated uniformly over the SDP. Data-driven corrections are applied to correct for known 
differences between data and simulation in the tracking, trigger and PID efficiencies, using 
identical methods to those described in Ref. 10 . The efficiency functions are obtained 


by fitting the corrected simulation with two-dimensional cubic splines. Figure [4] shows 
the histograms used to model the variation of the efficiency over the SDP for candidates 
triggered by (a) signal decays and (b) by the rest of the event. 


6.2 Background studies 

Among the background yields in the signal region, given in Table [lj the only sizeable 
components are due to combinatorial background and B° —» D ( '*" i0 tt + tt~ decays. The SDP 
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Figure 5: SDP distributions of the background contributions from (a) combinatorial background 
and (b) B° —> Zh*)°decays. 


distributions of these backgrounds are described by the histograms shown in Fig. [5j 

Combinatorial background contributes 22 % of candidates in the signal region. The 
shape of this contribution over the SDP is obtained from the high B mass sideband (5400- 
5900 MeV). As seen in Fig. [2j this sideband is dominated by combinatorial background, 
with a small contribution from B° — » decays. The SDP distribution of the 

B° — » ZT*)°7r + 7r - contribution in the sideband region is modelled as described below 
for the signal region and is subtracted from the sideband sample. To check that the 
SDP shape of the combinatorial background does not vary with the B candidate mass, 
a sample of doubly charged D°K ± tt ± candidates is investigated. This confirms that the 
SDP distribution from the sideband is a reliable model for combinatorial background in 
the signal region. 

The B° ->■ DW° 7r + 7r contribution accounts for 1.6% of the candidates in the signal 
region. Its SDP distribution is obtained from a simulated sample, with the B° —y D°7r + 7r _ 
and B° —* D*°7t + tt^ contributions combined and reweighted as described in Sec. [d] The 
B° —y decay mode is the dominant component of this background in the signal 

region. 

Due to resolution effects, B° —y D*(2010)~K + , ZT(2010) - — > D° decays are not 
entirely removed by the D* veto. Although this component corresponds to real B° —y 
D°n~K + decays, it is treated as incoherent (he. as background) since its interference with 
the remainder of the signal is negligible. Its SDP distribution is modelled as a Gaussian 
peak in m(D 0 7i~), with mean fixed to the known D*(2010)“ mass and width and yield 
determined in the fit to data. 
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Table 2: Signal contributions to the fit model, where parameters and uncertainties are taken 
from Ref. 


37 . The models are described in Sec. 5 


Resonance 

Spin 

DP axis 

Model 

Parameters 

K* (892)° 

1 

m 2 (K + TT~ ) 

RBW 

m 0 = 895.81 ± 0.19 MeV, r 0 = 47.4 ± 0.6 MeV 

K* (1410)° 

1 

m 2 (K + TT ~) 

RBW 

m 0 = 1414 ± 15 MeV, r 0 = 232 ± 21 MeV 

K* (1430)° 

0 

m 2 (K + n~) 

LASS 

Determined from data (see text) 

K* (1430)° 

2 

m 2 (K + TT ~) 

RBW 

m 0 = 1432.4 ± 1.3 MeV, r 0 = 109 ± 5 MeV 

Dq(2400)“ 

0 

m 2 (O 0 7r“) 

RBW 

Determined from data (see Table [3j 

0^(2460)“ 

2 

m 2 (D°TT~) 

RBW 


Nonresonant 

0 

m 2 (D°TT~) 

dabba 

Fixed (see text) 

Nonresonant 

1 

m 2 (D° 7r“) 

EFF 

Determined from data (see text) 


6.3 Amplitude model for B° —> D°K + ir decays 


The Dalitz plot fit is performed using the Laura++ [53] package. The likelihood function 
is given by 


N c 




y N t V t ),mJ(A'+7r )) 


(16) 


where the index i runs over N c candidates, k runs over the signal and background 
components and TV*, is the yield in each component. The probability density function for 
the signal component, V S i g , is given by Eq. (12) with the | A ^m 2 (D°7r~), m 2 (A' + 7r“)) | 2 
terms multiplied by the efficiency function described in Sec. |6.1 As it is possible for 


the minimisation procedure to find a local minimum of the negative log likelihood (NLL) 
profile, the fit is repeated many times with randomised initial values to ensure the global 
minimum is found. 

The nominal Dalitz plot fit. model for B° —>■ D°K + 7r~ decays is composed of several 
resonant and nonresonant. amplitudes. Only those amplitudes that provide significant 
contributions or that aid the fit stability are included. Unnatural spin-parity states are 
not considered as these do not decay to two pseudoscalars. The eight amplitudes included 
in this model are listed in Table [2] The width of the narrowest signal contribution to 
the Dalitz plot (~ 50MeV) is far larger than the mass resolution (~ 2.4 MeV); therefore, 
resolution effects are neglected. 

The real and imaginary parts of the complex coefficients Cj defined in Eq. (JTj) are free 
parameters of the fit except for the coefficient of the .02(2460)“ component, which is fixed 
to 1 as a reference. The phases and magnitudes of the complex coefficients, as well as fit 
fractions and interference fit fractions are derived from these free parameters. In addition, 
the masses and widths of the Dq( 2400)“ and .02(2460)“ resonances are determined from 
the fit to data and are reported in Table [3] The statistical uncertainties on all parameters 
of interest are calculated using large samples of simulated pseudoexperiments. These 
pseudoexperiments are also used to determine the correlations between the statistical 
uncertainties on the parameters, which are given in Appendix [A] The LASS parameters 
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Table 3: Masses and widths (MeV) determined in the fit to data, with statistical uncertainties 
only. 


Resonance 

Mass 

Width 

L>*(2400)- 

2360 ± 15 

255 ± 26 

£>*(2460)" 

2465.6 ± 1.8 

46.0 ±3.4 


Table 4: Complex coefficients and fit fractions determined from the Dalitz plot fit. Uncertainties 
are statistical only. Note that the fit fractions, magnitudes and phases are derived quantities. 


Resonance 

Fit fraction 


Isobar model coefficients 



(%) 

Real part 

Imaginary part 

Magnitude 

Phase 

K* (892)° 

37.4 ± 1.5 

-0.00 ±0.15 

-1.27 ±0.06 

1.27 ±0.06 

— 1.57 ± 0.11 

IC (1410)° 

0.7 ±0.3 

0.15 ±0.06 

-0.09 ±0.09 

0.18 ±0.07 

-0.54 ±0.21 

A* (1430)° 

5.1 ±2.0 

0.14 ±0.38 

0.45 ±0.15 

0.47 ±0.09 

1.27 ±0.95 

LASS nonresonant 

4.8 ±3.8 

-0.10 ±0.24 

0.44 ±0.14 

0.46 ±0.14 

1.79 ±0.65 

LASS total 

6.7 ± 2.7 





K*_ (1430)° 

7.4 ± 1.7 

-0.32 ±0.09 

-0.47 ±0.07 

0.57 ±0.05 

—2.16 ± 0.19 

£($(2400)- 

19.3 ±2.8 

-0.80 ±0.08 

-0.44 ±0.14 

0.91 ±0.07 

-2.64 ±0.15 

£>2(2460)- 

23.1 ±1.2 

1.00 

0.00 

1.00 

0.00 

Dir S-wave (dabba) 

6.6 ± 1.4 

-0.39 ±0.09 

0.36 ±0.17 

0.53 ±0.07 

2.40 ±0.27 

Dir P-wave (EFF) 

8.9 ± 1.6 

-0.62 ±0.06 

-0.03 ±0.06 

0.62 ±0.06 

-3.09 ±0.10 

Total fit fraction 

113.4 






are determined to be m 0 = 1450 ± 80 MeV, T 0 = 400 ± 230 MeV, a = 3.2 ± 1.8 GeV and 
r = 0.9 ±1.1 GeV, while the parameter of the EFF lineshape of the Dir P-wave nonresonant 
amplitude is determined to be a = 0.88 ± 0.10 GeV~ 2 . 

The values of the fit fractions and complex coefficients obtained from the fit are shown 
in Table [4j while the values of the interference fit fractions are given in Appendix [B] The 
sum of the fit fractions is seen to exceed unity, mostly due to interference within the Dir and 
Kir S-waves. Note that in Table |4| and all results for fit fractions and derived quantities, 
values are reported separately for both the Kq (1430)° and nonresonant components of the 
LASS lineshape, as well as their coherent sum. 

Projections of the data and the nominal fit model onto m{K + ir~), m(D°ir~) and 
m(D°K + ) are shown in Fig. [bj Zooms are provided in the regions of the main resonant 
contributions in Fig. [t} Projections onto the cosine of the hclicity angle of the D°ir~ and 
K + it~ systems are shown in Figs. [8] and [9j respectively. These projections all show good 
agreement between data and the fit model. 

Angular moments provide a useful method to investigate the helicity structure of the 
decays. The angular moments in m{D°ir~) and m(K + ir~), obtained as in Refs. 10,11 , 
are shown in Figs. 10 and 11, respectively. The contributions due to the K* (892)° and 
-D?; (2460) - resonances are seen as peaks in moments up to order 2 and order 4, respectively, 
as expected for spin-1 and spin-2 resonances. Reflections make the interpretation of 
moments at higher masses more difficult. However, the good agreement seen between data 
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and the fit model provides further support for the fit model being a good description of 
the data. 

The quality of the fit is evaluated by determining a two-dimensional x' 2 value, comparing 
the data and fit model in f44 bins across the SDP that are defined adaptively to ensure 
they are approximately equally populated. The pull in each of these bins, defined as the 
difference between the data and the fit model divided by the uncertainty, is shown in 
Fig. 12 The effective number of degrees of freedom of the \ 2 is between iVbins 


iV, 


pars 


and iVbins— 1, where N pars is the number of free parameters in the fit, yielding a reduced \ 2 
in the range 0.99 - 1.22. Additional unbinned tests of fit quality 54 also show acceptable 
agreement between the data and the fit model. 


7 Systematic uncertainties 


Systematic uncertainties are considered from sources that can be divided into two categories, 
experimental and model uncertainties. Experimental uncertainties arise due to the signal 
and background yields in the signal region, the distributions of the background components 
across the SDP, the variation of efficiency across the SDP, and possible bias due to the 
fitting procedure. Model uncertainties are considered due to the fixed parameters in the 
signal model, the addition or removal of marginal amplitudes, and the choice of models for 
the Kit S-wave and the Dtt S- and P-waves. The uncertainties due to all of these sources 
are combined in quadrature. 

The signal and background yields are determined from the fit to the B candidate 
invariant mass distribution. Systematic uncertainties are considered on the total yields 
due to both statistical and systematic variations of these yields, as evaluated in Ref. 15 


The yields in the signal region are varied appropriately and the effects of these variations 
on the parameters of interest are assigned as systematic uncertainties. 

The histograms describing the background distributions are each varied within statistical 
uncertainties to establish the uncertainty due to the imperfect knowledge of the distribution 
of these backgrounds across the SDP. The histograms describing the efficiency variation 
across the SDP are also varied within their uncertainties prior to the spline fit. The 
efficiency histograms are also varied with local correlations to account for the effects of 
any unknown correlations between neighbouring bins. For each parameter the larger 
uncertainty is considered from the correlated and uncorrelated variations. In addition, the 
binning scheme of the control sample used to evaluate PID performance is varied. 

An ensemble of pseudoexperiments is used to establish systematic uncertainties related 
to any potential fit bias. Pseudoexperiments are generated using the parameters determined 
in the fit to data, and mean fitted values of the parameters are determined from Gaussian 
fits to the distributions of results obtained in the ensemble. The differences between the 
generated and mean fitted values are found to be small. Uncertainties are determined as 
the sum in quadrature of the difference between the generated and mean fitted values, and 
the uncertainty on the mean fitted value from the Gaussian fit. 

Systematic uncertainties due to fixed parameters in the fit model are determined by 
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Data 

Kji S-wave (LASS) 
Dji S-wave (EFF) 


Total fit 


1 K*( 892)° - K*(1410)° 


- K*(1430)° 

- Dn P-wave (EFF) 


D * (2400) - D* (2460) 

D* Background 


Figure 6: Projections of the data and amplitude fit results onto (a) m(D°TT~), (c) vi{K + tt~) 
and (e) m(D°K + ), with the same projections shown in (b), (d) and (f) with a logarithmic y -axis 
scale. Components are described in the legend. 
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Figure 7: Projections of the data and amplitude fit results onto (a) m(D°n ) in the L>2(2460) 
region and (b) the low m(A + 7r _ ) region. Components are as shown in Fig.pl 




Figure 8: Projections of the data and amplitude fit results onto cos0(Z?°7r _ ) in the mass ranges 
(a) 2.04 < m(D°TT~) < 2.35 GeV and (b) 2.35 < m(D°- k~) < 2.55 GeV. Components are as 
shown in Fig. [6j 


repeating the fit with these parameters varied. The fixed masses and widths are varied 
within their uncertainties as shown in Table[2j while the Blatt-Weisskopf barrier radii,_jBW 7 


are varied independently for D°7r~ and K + n~ resonances between 3 and 5 GeV -1 1101. 

The K*{ 1410)° amplitude does not contribute significantly to the nominal fit model 
so the fit is repeated with this component removed and the change in results is assigned 
as the associated systematic uncertainty. In addition, uncertainties are assigned from 
the changes in the results when an amplitude due to the A'*(1680)° resonance or from a 
virtual A* (5325) 


component, as described in Ref. 10 


A Dalitz plot analysis of B° 
at m(D°K~) ~ 2.86GeV |9, 
2.76 GeV in B° ->■ D° tt+tt- 


D°K~ 


7 T 1 


10 


is included in the model, 
observed both spin-1 and spin-3 resonances 


spin-1 resonance at similar m(D +/ n ) values 


A spin-3 resonance has also been seen at m(D° ti~) ~ 
decays 112], while an analysis of B + —> D~K + ti + shows a 
Uncertainties are assigned corresponding 
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Figure 9: Projections of the data and amplitude fit results onto cos 0{K + it~) in the mass ranges 
(a) m(K + 7r~) < 0.8GeV, (b) 0.8 < m{K + ir~) < 1.0GeV, (c) 1.0 < m(K + n~) < 1.3GeV and 
(d) 1.4 < m(K + ir~) < 1.5 GeV. Components are as shown in Fig. [gJ 


to the changes in results when either a spin-1 or spin-3 state at m(D°7T ~) ~ 2.76 GeV is 
added to the nominal model. 

The LASS model used to describe the Ktt S-wave is replaced with a Flatte shape [55] , 
which accounts for the Krf threshold near the A'q( 1430)° mass, and a resonant term with 
a modified mass dependent width 52 for the k (also known as A"o(800)) state. This 
model gives an NLL that is worse by 4.4 units. A model-independent description of the 
Ktt S-wave is also used to fit the data, with the uncertainty on each parameter taken 
as the larger of the two differences. This alternative improves the NLL by 8.8 units but 
introduces an additional four free parameters into the fit. The dabba lineshape used to 
describe the D°tt~ S-wave is replaced with an EFF lineshape, while a power-law model is 
introduced as an alternative to the EFF lineshape that describes the D°tt~ P-wave. These 
alterations worsen the NLL by 3.0 units and improve it by 0.8 units respectively. 

The total systematic uncertainties due to experimental and model effects for fit fractions 
and complex coefficients are given together with the results in Sec. [8] The contributions 
from different sources to the systematic uncertainties on the fit fractions, masses and 
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Figure 10: Background-subtracted and efficiency-corrected Legendre moments up to order 7 
calculated as a function of m(D 0/ K~) for data (black data points) and the fit result (solid blue 
curve). 
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Figure 11: Background-subtracted and efficiency-corrected Legendre moments up to order 7 
calculated as a function of m(K + ir~) for data (black data points) and the fit result (solid blue 
curve). 
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Figure 12: Differences between the data SDP distribution and the fit model across the SDP, in 
terms of the per-bin pull. 


Table 5: Experimental systematic uncertainties on the fit fractions (%) and masses and widths 
(MeV). Uncertainties given on the central values are statistical only. 


Resonance 

Central value 

S/B fraction 

Efficiency 

Background SDP 

Fit bias 

Total 

K* (892)° 

37.4 ± 1.5 

0.60 

0.83 

0.50 

0.31 

1.17 

AT*(1410)° 

0.7 ± 0.3 

0.06 

0.39 

0.69 

0.05 

0.80 

A'q(1430)° 

5.1 ±2.0 

0.28 

1.48 

1.85 

0.33 

2.41 

LASS nonresonant 

4.8 ±3.8 

0.51 

2.25 

2.86 

0.86 

3.77 

LASS total 

6.7 ± 2.7 

0.26 

1.86 

1.60 

1.02 

2.67 

AT 2 *(1430)° 

7.4 ± 1.7 

0.23 

0.72 

0.53 

0.54 

1.07 

£>2(2400)- 

19.3 ±2.8 

0.21 

1.39 

1.43 

0.40 

2.04 

£>2(2460)“ 

23.1 ± 1.2 

0.70 

0.70 

0.49 

0.15 

1.11 

£>7t S-wave (dabba) 

6.6 ± 1.4 

0.03 

0.81 

0.59 

0.57 

1.15 

£>7t P-wave (EFF) 

8.9 ± 1.6 

0.86 

1.91 

0.52 

0.38 

2.19 

m (£>2(2400)-) 

2360 ± 15 

4.6 

8.1 

7.0 

3.7 

12.2 

m (£>2(2460)-) 

2465.6 ± 1.8 

0.01 

0.37 

0.22 

0.29 

0.51 

T (£>2(2400)-) 

255 ± 26 

2.8 

13.1 

13.9 

4.8 

19.9 

r (£>2(2460)-) 

46.0 ±3.4 

0.5 

0.9 

0.9 

0.5 

1.4 


widths are given in Tables [5] and [6j The dominant experimental systematic uncertainty 
on most parameters is due to either the efficiency variation or background distributions 
across the SDP, while the model uncertainties are generally dominated by the effects of 
alternative descriptions of the Kit S-wave and the Dn S- and P-waves. For the parameters 
of broad components, the model uncertainties dominate; however for other parameters the 
statistical, experimental and model uncertainties are all of a similar magnitude. 

The stability of the fit results is confirmed by several cross-checks. The data sample is 
divided into subsamples according to the flavour of the B candidate, the trigger decision, 
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Table 6: Model uncertainties on the fit fractions (%) and masses and widths (MeV). Uncertainties 
given on the central values are statistical only. 


Resonance 

Central value 

Fixed parameters 

Add/remove 

Alternative models 

Total 

K*( 892)° 

37.4 ± 1.5 

0.75 

1.14 

1.09 

1.74 

£T*(1410)° 

0.7 ±0.3 

0.18 

0.70 

0.22 

0.76 

K* (1430)° 

5.1 ±2.0 

0.79 

3.30 

0.23 

3.40 

LASS nonresonant 

4.8 ±3.8 

1.10 

3.99 

5.20 

6.65 

LASS total 

6.7 ±2.7 

0.53 

1.42 

5.21 

5.43 

K* (1430)° 

7.4 ± 1.7 

0.36 

1.87 

0.56 

1.98 

Ug (2400)“ 

19.3 ±2.8 

0.55 

1.95 

7.11 

7.40 

£>2(2460)- 

23.1 ± 1.2 

0.18 

0.73 

0.99 

1.24 

Dtt S-wave (dabba) 

6.6 ± 1.4 

0.27 

1.40 

3.46 

3.74 

Dtt P-wave (EFF) 

8.9 ± 1.6 

0.31 

1.99 

2.15 

2.95 

m (£>2(2400)-) 

2360 ± 15 

6.1 

9.3 

25.6 

27.9 

??z (£>3(2460)-) 

2465.6 ± 1.8 

0.09 

1.05 

0.48 

1.15 

r (£>2(2400)-) 

255 ± 26 

4.0 

18.0 

43.5 

47.2 

r (£>2(2460)-) 

46.0 ±3.4 

1.4 

0.5 

2.4 

2.9 


the polarity of the magnet and the year of data taking. Each subsample is fitted separately 
and all are seen to be consistent with the nominal fit results. The fit is also repeated 
with additional resonance components to test the fit model. Resonances up to spin-3 were 
considered for all three pairs of daughters with the mass and width allowed to vary in the 
fit. No additional significant contributions were observed. 


8 Results and summary 


The results for the complex coefficients are reported in Table [7J with results for the 
fit fractions given in Table [8j The results for the interference fit fractions are given 
in Appendix [b] Upper limits are determined on the fit fractions of the AT* (1410)° and 
D^(2760)~ components, which do not give significant contributions, using the method 
described in Ref. [l0|. 

The fit fractions for resonant contributions are converted into quasi-two-body product 
branching fractions by multiplying by B(B° —>• D 0 K + tt~) = (9.2 ± 0.6 ± 0.7 ± 0.6) x 10 -5 , 
which is obtained from the ratio of branching fractions B(B° —>■ D°K + 7r~)/B(B 0 
D° tt h 


—>■ 


7T 


15 


multiplied by the latest result for B(B° —» H°7r + 7r _ ) |12], accounting for 
the different £>*(2010)- veto windows used in the analyses and for the 3.7% residual 
background due to B° — y D*(2010)~K + decays. The statistical correlation between the 
results of Ref. 15 and this analysis is found to be negligible. The results for the product 
branching fractions are shown in Table [9| For the K + n~ resonances, where the branching 
fractions to Kn are known, the product branching fractions are converted to the B 

The results for B(B° — » D°K*( 892)°) and 
are consistent with, and more precise 


10 


decay branching fractions, shown in Table 
B(B° -± £>*(2460)-/!+) x £(£>*(2460)- -± £>°tT 
than, previous measurements [14,16 \m- 
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Table 7: Results for the complex amplitudes and their uncertainties presented (top) in terms of 
real and imaginary parts and (bottom) in terms and magnitudes and phases. The three quoted 
errors are statistical, experimental systematic and model uncertainties, respectively. 


Resonance 


Real part 


Imaginary part 


K*( 892)° 
Jl*(1410)° 

Jig (1430)° 

LASS nonresonant 
IP*(1430)° 

D*( 2400)- 
D*(2460)“ 


-0.00 ±0.15 ±0.24 ±0.34 
0.15 ±0.06 ±0.04 ±0.09 
0.14 ±0.38 ±0.48 ±0.38 
-0.10 ±0.24 ±0.16 ±0.42 
-0.32 ±0.09 ±0.15 ±0.23 
-0.80 ±0.08 ±0.07 ±0.22 
1.00 


-1.27 ±0.06 ±0.03 ±0.06 
-0.09 ±0.09 ±0.18 ±0.18 
0.45 ±0.15 ±0.37 ±0.17 
0.44 ±0.14 ±0.17 ±0.23 
-0.47 ±0.07 ±0.14 ±0.15 
-0.44 ±0.14 ±0.12 ±0.18 
0.00 


Dn S-wave 
Dn P-wave 

(dabba) 

(EFF) 

-0.39 ±0.09 ±0.09 ±0.14 

-0.62 ±0.06 ±0.03 ±0.11 

0.36 ±0.17 ±0.14 ±0.23 

-0.03 ±0.06 ±0.05 ±0.10 


Resonance 


Magnitude 

Phase 


K* (892)° 

K* (1410)° 

K* (1430)° 

LASS nonresonant 
iP*(1430)° 
Hq(2400) _ 
±>^(2460)- 


1.27 ± 0.06 ± 0.03 ± 0.05 -1.57 ± 0.11 ± 0.16 ± 0.27 

0.18 ± 0.07 ± 0.10 ± 0.11 -0.54 ± 0.21 ± 0.55 ± 1.04 

0.47 ± 0.09 ± 0.10 ± 0.14 1.27 ± 0.95 ± 1.04 ± 0.81 

0.46 ± 0.14 ± 0.16 ± 0.29 1.79 ± 0.65 ± 0.35 ± 0.69 

0.57 ± 0.05 ± 0.04 ± 0.08 -2.16 ± 0.19 ± 0.43 ± 0.43 

0.91 ± 0.07 ± 0.06 ± 0.17 -2.64 ± 0.15 ± 0.14 ± 0.23 

1.00 0.00 


Dn S-wave (dabba) 
Dn P-wave (EFF) 


0.53 ± 0.07 ± 0.04 ± 0.14 2.40 ± 0.27 ± 0.24 ± 0.44 
0.62 ± 0.06 ± 0.04 ± 0.11 -3.09 ± 0.10 ± 0.07 ± 0.17 


Table 8: Results for the fit fractions and their uncertainties (%). The three quoted errors are 
statistical, experimental systematic and model uncertainties, respectively. Upper limits are given 
at 90 % (95 %) confidence level. 


Resonance 

Fit fraction 

LIpper limit 

K* (892)° 

37.4 ± 1.5 ± 1.2 ± 1.7 


K* (1410)° 

0.7 ±0.3 ±0.8 ±0.8 

< 3.2 (3.7) 

A'g(1430)° 

5.1 ±2.0 ±2.4 ±3.4 


LASS nonresonant 

4.8 ±3.8 ±3.8 ±6.7 


LASS total 

6.7 ±2.7 ±2.7 ±5.4 


JP|(1430)° 

7.4 ± 1.7 ± 1.1 ±2.0 


D* (2400)- 

19.3 ±2.8 ±2.0 ±7.4 


±>*(2460)- 

±>*(2760)“ 

23.1 ±1.2 ±1.1 ±1.2 

< 1.0 (1.1) 

Dn S-wave (dabba) 

6.6 ±1.4 ±1.2 ±3.7 


Dir P-wave (EFF) 

8.9 ±1.6 ±2.2 ±3.0 
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Table 9: Results for the product branching fractions. The four quoted errors are statistical, 
experimental systematic, model and PDG uncertainties, respectively. Upper limits are given at 
90 % (95 %) confidence level. 


Resonance 

Product branching fraction (10 5 ) 

Upper limit (10 5 ) 

K*( 892)° 

3.42 ± 0.13 ± 0.10 ± 0.16 ± 

0.40 


K* (1410)° 

0.07 ± 0.03 ± 0.08 ± 0.07 ± 

0.01 

< 0.29 (0.34) 

o * 

CO 

o 

o 

0.47 ± 0.18 ± 0.22 ± 0.31 ± 

0.05 


LASS nonresonant 

0.44 ±0.34 ±0.34 ±0.61 ± 

0.05 


LASS total 

0.61 ± 0.25 ± 0.25 ± 0.49 ± 

0.07 


K* (1430)° 

0.68 ± 0.15 ± 0.10 ± 0.18 ± 

0.08 


£>*(2400)- 

1.77 ± 0.26 ± 0.19 ± 0.67 ± 

0.20 


£>2(2460)“ 

£>3(2760)- 

2.12 ± 0.10 ± 0.11 ± 0.11 ± 

0.25 

< 0.10 (0.11) 

£>7r S-wave (dabba) 

0.60 ± 0.13 ± 0.11 ± 0.34 ± 

0.07 


Dtt P-wave (EFF) 

0.81 ± 0.15 ± 0.20 ± 0.27 ± 

0.09 



Table 10: Results for the branching fractions. The four quoted errors are statistical, experimental 
systematic, model and PDG uncertainties, respectively. Upper limits are given at 90 % (95 %) 
confidence level. 


Resonance 

Branching fraction (10 E 

’) Upper limit (10 5 ) 

K* (892)° 

5.13 

± 

0.20 

± 

0.15 

± 

0.24 

± 

0.60 


K* (1410)° 

1.59 

± 

0.68 

± 

1.81 

± 

1.59 

± 

0.36 

< 6.7 (7.8) 

K* (1430)° 

0.71 

± 

0.27 

± 

0.33 

± 

0.47 

± 

0.08 


LASS non resonant 

0.66 

± 

0.51 

± 

0.51 

± 

0.92 

± 

0.08 


LASS total 

0.92 

± 

0.38 

± 

0.38 

± 

0.74 

± 

0.11 


A* (1430)° 

2.04 

± 

0.45 

± 

0.30 

± 

0.54 

± 

0.25 



The masses and widths of the £>q( 2400) and £>2(2460) states are found to be 


m(D*( 2400)") 
r(£>*(2400)-) 
m(D*( 2460)“) 
r(L>;(2460)“) 


(2360 ± 15 ± 12 ± 28) MeV 
(255 ± 26 ± 20 ± 47) MeV 
(2465.6 ± 1.8 ± 0.5 ± 1.2) MeV 
(46.0 ± 3.4 ± 1.4 ±2.9) MeV, 


where the three uncertainties are statistical, experimental systematic and model systematic, 
respectively. These are consistent with, though less precise than, recent results from a DP 
analysis of B° —> D °decays 12 . They also show good agreement with, and have 


13,37 


similar precision to, earlier measurements of these quantities 

In summary, the first amplitude analysis of B° —> D°K + ir~ decays has been presented, 
using a data sample corresponding to 3.0 RU 1 of pp collision data collected by the LHCb 
experiment. A good description of the data is obtained with a model containing contri- 
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butions from intermediate K*( 892)°, A'* (1410)°, A"2(1430) 0 and Z^2(2460) _ resonances, 
with additional components to describe broad structures in the Kit S-wave and the Dit S- 
and P-waves. The masses and widths of the .Dq( 2400) _ and Ah; (2460) ~ resonances are 
measured, as are the complex amplitudes and fit fractions for all components included 
in the amplitude model. The results can be used in conjunction with an analysis of 
B° —> DK + it~ decays, where the neutral D meson is reconstructed in final states such as 
K + K ~, to measure the CKM unitarity triangle parameter 7 [l, 21. 
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Appendices 

A Correlation matrices 

The statistical correlations between the complex coefficients, masses and widths determined 
from the fit to data are given in terms of real and imaginary parts, and in terms of 
magnitudes and phases in Tables 0 and [12} The correlations are determined using the 
same sample of simulated pseudoexperiments used to calculate the statistical uncertainties 
on fit parameters. 


B Results for interference fit fractions 


The central values of the interference fit fractions are given in Table 13 The statist! 


cal, experimental systematic and model uncertainties on these quantities are given in 
Tables [T4l 


15 and 16 
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Table 13: Interference fit fractions (%) from the nominal DP fit. The amplitudes are all pairwise 
products involving: (Ao) A"*(892)°, (Ai) A'*(1410)°, (A 2 ) A'q(1430)°, (A 3 ) LASS nonresonant, 
(A 4 ) Ji2(1430)°, (A 5 ) Dq( 2400) _ , (Ag) Ahj(2460) _ , (A 7 ) Dir S-wave (dabba), (Ag) Dir P-wave 
(EFF). The diagonal elements are the same as the conventional fit fractions. 



Ao 

At 

A2 

A3 

A 4 

A5 

Aq 

a 7 

As 

Aq 

37.4 

1.8 

- 0.0 

0.0 

0.0 

-2.4 

- 0.2 

- 2.6 

-5.0 

A 1 


0.7 

- 0.0 

0.0 

0.0 

- 1.1 

-0.3 

- 0.6 

- 1.1 

A-2 



5.1 

-3.2 

- 0.0 

2.9 

1.5 

1.6 

1.9 

^3 




4.8 

0.0 

-8.4 

0.1 

2.1 

0.2 

A 4 





7.4 

-0.3 

- 1.0 

-1.5 

-0.4 

A5 






19.3 

- 0.0 

2.8 

0.0 

Aq 







23.1 

- 0.0 

0.0 

A 7 








6.6 

- 0.0 

^8 









8.9 


Table 14: Statistical uncertainties on the interference fit fractions (%). The amplitudes are 
all pairwise products involving: (Ao) A"*(892)°, (Ai) /l*(1410)°, (A 2 ) Ad(1430)°, (A 3 ) LASS 
nonresonant, (A 4 ) Jl|(1430)°, (A 5 ) A>d(2400) _ , (Ag) ZT;(2460) - , (A 7 ) Dir S-wave (dabba), (Ag) 
Dir P-wave (EFF). The diagonal elements are the same as the conventional fit fractions. 



Ao 

A 1 

A 2 

A 3 

7U 

A 5 

Ao 

A 7 

^8 

Ao 

1.5 

0.7 

0.0 

0.0 

0.0 

0.9 

0.4 

0.4 

0.5 

A 1 


0.3 

0.0 

0.0 

0.0 

0.7 

0.3 

0.4 

0.5 

A 2 



2.0 

2.6 

0.0 

1.1 

0.5 

1.5 

0.5 

A 3 




3.8 

0.0 

3.2 

0.4 

3.4 

0.8 

a 4 





1.7 

0.6 

0.1 

0.6 

0.5 

A 5 






2.8 

0.0 

3.1 

0.0 

Ao 







1.2 

0.0 

0.0 

A 7 








1.4 

0.0 

As 









1.6 
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Table 15: Experimental systematic uncertainties on the interference fit fractions (%). The 
amplitudes are all pairwise products involving: (Ao) 7i*(892)°, (A\) A'*(1410)°, (A 2 ) 4vq(1430)°, 
(A 3 ) LASS nonresonant, (A 4 ) Ji2(1430)°, (A 5 ) 44q( 2400)“, (Ag) 442 ( 2460 )“, (A 7 ) Dir S-wave 
(dabba), (As) Dir P-wave (EFF). The diagonal elements are the same as the conventional fit 
fractions. 



Ao 

A± 

A2 

^3 

A4 

A5 

Aq 

a 7 

7^8 

Ao 

1.2 

0.5 

0.0 

0.0 

0.0 

2.3 

0.2 

0.4 

0.5 

Ai 


0.8 

0.0 

0.0 

0.0 

0.4 

0.8 

0.7 

0.9 

A2 



2.4 

3.1 

0.0 

1.4 

2.0 

2.6 

1.1 

A3 




3.8 

0.0 

3.6 

0.2 

2.6 

0.3 

a 4 





1.1 

2.2 

0.2 

0.5 

1.0 

A5 






2.0 

0.0 

4.0 

0.0 

Aq 







1.1 

0.0 

0.0 

a 7 








1.2 

0.0 

As 
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Table 16: Model uncertainties on the interference fit fractions (%). The amplitudes are all pairwise 
products involving: (Ao) A'*(892)°, (Ai) 4i*(1410)°, (A 2 ) A'q(1430)°, (A 3 ) LASS nonresonant, 
(A 4 ) Ji2(1430)°, (A 5 ) L>q( 2400) - , (Ag) 442 ( 2460 )“, (A 7 ) Dir S-wave (dabba), (Ag) Dir P-wave 
(EFF). The diagonal elements are the same as the conventional fit fractions. 



Ao 

Ai 

A2 

7h 

A4 

A 5 

Ao 

a 7 

A% 

Ao 

1.7 

0.8 

0.0 

0.0 

0.0 

2.6 

0.7 

0.8 

1.4 

a 4 


0.8 

0.0 

0.0 

0.0 

0.8 

0.6 

0.8 

1.1 

A2 



3.4 

2.9 

0.0 

5.8 

1.7 

3.6 

2.7 

A3 




6.6 

0.0 

8.3 

0.3 

3.5 

0.7 

a 4 





2.0 

1.7 

0.2 

1.2 

1.0 

A5 






7.4 

0.0 

4.5 

0.0 

A% 







1.2 

0.0 

0.0 

a 7 








3.7 

0.0 

7^8 









2.9 
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